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Abstract. The physical interpretation of spectropolarimetric observations 
of lines of neutral helium, such as those of the 10830 A multiplet, represents an 
excellent opportunity for investigating the magnetism of plasma structures in the 
solar chromosphere. Here we present a powerful forward modeling and inversion 
code that permits either to calculate the emergent intensity and polarization for 
any given magnetic field vector or to infer the dynamical and magnetic properties 
from the observed Stokes profiles. This diagnostic tool is based on the quantum 
theory of spectral line polarization, which self-consistently accounts for the Hanle 
and Zeeman effects in the most general case of the incomplete Paschen-Back 
effect regime. We also take into account radiative transfer effects. An efficient 
numerical scheme based on global optimization methods has been applied. Our 
Stokes inversion code permits a fast and reliable determination of the global 
minimum. 



1. Introduction 

The quality of spectropolarimetric observations of solar chromospheric plas- 
mas is steadily increasing. However, in order to obtain reliable empirical in- 
formation on the strength and geometry of the magnetic field vector we also 
need to develop and apply suitable diagnostic tools within the framework of 
th e quantum theory of spectral line pola rization (e.g., the recent monograph 
by lLandi Degl'Innocenti &: Landolfi 2004). In particular, the development of a 



forward modeling and inversion code of Stokes profiles is of great importance 
because it would facilitate the investigation of a variety of problems of great 
astrophysical interest. For such a diagnostic tool to be reliable, it has to include 
all the important physical mechanisms that in duce spectral line pol arization in 
stellar atmospheres (e.g., the recent review bv lTruiillo Buenol I2006T ) . Not only 
the Zeeman effect, but also the atomic level polarization induced by anisotropic 
radiation pumping and its modification by the Hanle effect have to be taken into 
account. For a remarkable example of a pioneering work in this field see Landi 
Degl'Innocenti (1982). 

One of the main motivations of our investigation has been to develop a 
robust but user-friendly computer program for facilitating the analysis of the 
polarization signals induced by the above-mentioned physical mechanisms in the 
spectral lines of the He I 10830 A multiplet, which are sens itive to both the Zee- 



man e ffect and to the presence of atomic level polarization (|Trujillo Bueno et al 
120021 ; iTruiillo Bueno Asensio Ramosl 120071 ). Spectro-polarimetric observa- 
tions of this multiplet have been carried out in a variety of plasma struc- 
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tures of the solar chromosphere and coron a, such as sunpots (Har vey Sz Halll 
19711: iRuedi et alJll995l:ICenteno et al.l l2006f). coronal filaments (iLin et alJll99£ " 



Trui illo Bueno et al.ll2002l ). prominences Truiillo Bu eno et al.ll200 2; Mcre nda et al 
2006), emerging magn etic flux regions (Solanki et al. 20031 : lLagg et al.l 12004! ) 



chromospheric spicules (jTruiillo Bueno et al. 20051 ). active region filaments (Martinez 
Pillet et al. 2007; in preparation) and flaring regions (Sasso et al. on p. 467 ff of 
these proceedings). 



2. The Forward Modeling Option 

The adopted atomic model includes the following five terms of the triplet system 
of He I: 2s 3 S, 2p 3 P, 3s 3 S, 3s 3 P and 3d 3 D. It has been concluded that this 11 
J-levels model with 6 radiative transitions between the terms is enough for a 
satisfactory modeling of the polarization properties of the spectral line s of the 
He I D 3 multiplet at 5876 A (|Bommierlll977l : lLandi Degl'Innocentilll982l ). and it 



should be also a reasonable approximation for those of the He 1 10830 A multiplet 
which result from transitions between the terms 2s 3 S and 2p 3 P. 

It is well-known tha t the He I atom can be c orrectly described by the L- 
S coupling scheme (e.g., ICondon Shortlevl Il935l ). The different J-levels are 



grouped in terms with well defined values of the electronic angular momen- 
tum L and the spin S. The energy separation between the J-levels inside each 
term is very small in comparison with the energy difference between different 
terms. Therefore, in addition to the atomic polarization of each J-level (popu- 
lation imbalances and coherences between its magnetic substates), it turns out 
to be fundamental to allow for coherences between different J-levels pertaining 
to the same term. However, coherences between different J-levels pertaining 
to different terms can be safely neglected. In conclusion, we can describe the 
atom via the formalism of the multi- term atom (see Sections 7.5 and 7.6 of 
Landi Deri'Innocenti fe Landolfi! l20bll ). 



The atomic level polarization is quantified by the multipole components, 
Pq(J, J ), of the atomic density matrix ( Landi Degl'Innocenti Sz Landolfi 2004h . 



The number of real quantities required to describe the excitation state of our 
multi-term model atom is 405. The helium atoms are illuminated by the (given) 
aniso tropic radiation field coming from the underlying solar photosphere (jPiercd 
2000), whose properties are quantified by two spheri cal tensors: Jq and Jq (see 



Eq. 5.164 in lLandi Degl'Innocenti &: Landolfi! 12004 ) . They describe the mean 



intensity and the "degree of anisotropy" of the radiation field, respectively. We 
also consider that a magnetic field B is present, with strength B, inclination 0b 
with respect to the local vertical and azimuth xb- The density matrix elements 
ar e obtained by solving the statistical e quilibrium equations (see Section 7.6a 
of Landi Degl'Innocenti Landolfi 2004! ) . which take into account the effect of 



the radiative transitions in the presence of the Zeeman splittings produced by 
the assumed magnetic field vector. 

Once the elements of the atomic density matrix are known, the coefficients 
of the emission vector and of the propagation matrix o f the Stokes- vector transfer 



equa tion can be directly calculated (see Section 7.6b of lLandi Degl'Innocenti Landolfi 



20041 ) . The next step is to compute the emergent Stokes profiles by solving the 



Stokes-vector transfer equation. To this aim, our code includes several options. 
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The simplest one assumes that the emergent Stokes profiles are simply propor- 
tional to the corresponding emissivity. This optically-thin case is representative 
of the conditions found in some prominences. The second option assumes that 
the He I atoms are located in a slab of constant physical properties character- 
ized by a non-negligible optical depth At and that magneto-optical effects are 
negligible (see Trujillo Bueno et al. 2005 for the analytical expressions of the 
emergent Stokes profiles). Finally, the most gen eral radiative transfe r option is 
accomplished by using the DELOPAR method (|Truiillo Buendl2003h . 

Following the previously described approach, we have developed an effi- 
cient computer program that we have combined with an easy-to-use interactive 
front-end. It allows to interactively select the physical effects to be taken into 
account, thus helping to achieve a fast and reliable investigation of the influ- 
ence of various physic al parameters on the emergent Stokes parameters (for 
a first application, see iTruiillo Bueno Asensio Ramosl 20071 ). For instance, 
our code allows to include or discard the effect of the magneto-optical terms 
and/or stimulated emission processes. Furthermore, it is possible to carry out 
the forward modeling including or neglecting the presence of atomic polariza- 
tion. When atomic polarization is not taken into account, the problem reduces 
to the simplified case of the Zeeman effect, but with the positions and strengths 
of the cr and tt components calculated within the framework of the incomplete 
Paschen-Back effect theory. This is important for the He I 10830 A multiplet 
because the li near Zeeman effect theory underestimates the inferred magnetic 
field strength ( Socas-Navarro et al. 20041 ) . In the absence of atomic level po- 



larization, the emergent spectral line radiation is polarized due to the non-zero 
magnetic splitting produced by the presence of the magnetic field. For com- 
pleteness, the Milne-Eddington solution of the radiative transfer equation can 
also be used in this case. By default, the total Hamiltonian (including the fine 
structure and magnetic contributions) is diagonalized numerically. However, our 
code also permits to calculate the emergent Stokes profiles assuming the linear 
Zeeman regime, which is useful in order to compare with those calculated within 
the framework of incomplete Paschen-Back theory. 



3. The Stokes Inversion Option 

The previously described forward modeling code has been also combined with a 
robust inversion scheme. The inversion algorithm is based on the minimization 
of a merit function that describes how well our model reproduces the observed 
Stokes profiles. This merit function is chosen to be the standard x 2 ~f unc tion 
(least square method). The minimization algorithm tries to find the value of the 
model's parameters that lead to synthetic Stokes profiles in closest agreement 
with the observed ones (i.e., to those that imply the smallest value of the merit 
function). The model's parameters are the magnetic field vector B, the thermal 
velocity vth, the macroscopic velocity v mac and the damping of the line a. If 
radiative transfer effects are taken into account, this set is augmented with the 
optical depth of the slab (At). 

The standard Levenberg-Marquardt (LM) procedure for the minimization 
of the x 2 function has been used. The LM method is one of the fastest and 
simplest available methods when the initial estimation of the parameters is close 
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Figure 1. This figure illustrates the philosophy of the DIRECT method for 
searching for the region where the global minimum is located. In this case, 
we illustrate a case in two dimensions. After the evaluation of the function at 
some selected points inside each region, the DIRECT algorithm decides, using 
the Lipschitz condition, which rectangles can be further subdivided because 
they are potentially interesting. The method rapidly locates the region where 
the minimum is located. 



to the correct minimum. On the contrary, like with the majority of the standard 
numerical methods for function minimization, its main drawback is that it can 
easily get trapped in a local minimum of the \ 2 function. A straightforward but 
time consuming strategy to overcome this difficulty is to restart the minimization 
process with different values of the initial parameters. 

In our case, this difficulty has been overcome by using global optimization 
methods (GOM) that provide a good starting point. The majority of this type 
of methods are based on stochastic optimization techniques. The key idea is to 
efficiently sample the whole space of parameters to find the global minimum of a 
given function. One of the most promising methods is genetic optimization. In 
spite of the lack of a c onvergence theorem, they perform quite well in practice 
(e.g., iLagg et aill2004l ). In our case, we propose to apply instead the DIRECT 



deterministic method that relies o n a strong mathem atical basis. The name 
stands for "Dividing RECTangles" ( Jones et al. 1993) and the idea is to recur- 



sively sample parts of the space of parameters, locating in each iteration the part 
of the space where the global minimum is potentially located. The decision algo- 
rithm is based on the assumption that the function is Lipschitz continuous. For 
clarifying purposes, consider the case of a one-dimensional function / : M — > R. 
It is said to be Lipschitz continuous with constant a in an interval M = [a, b] if: 

\f(x) - f{x')\ < a\x -x'\, Vx,x'GM. (1) 

In particular, a consequence of the previous equation is that the two following 
equations have to be fulfilled Vx G M: 

Six) > f(a)-a(x-a) (2) 
/(*) > f(b) + a(x-b). (3) 

These two straight lines form a V-shape below the function f(x) whose intersec- 
tion provides a first estimation x\ of the minimum of the function. Repetition 
of this procedure in the two new subintervals [a, x%] and [xi,b] leads to new es- 
timations Xi that converge to the global minimum. Severe problems arise when 
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Step 


Method 


Free parameters 


Stokes profiles 


1 


DIRECT 


vth, v mac , At, a 


/ 


2 


LM 


Vth, "mac, At, a 




3 


DIRECT 


B, 9 B , Xb 


i, Q, u, v 


4 


LM 


B, 9b, Xb 


i, Q, u, v 



Table 1. Scheme applied for the inversion of the Stokes profiles. The two 
DIRECT cycles are used for obtaining an initial value of the parameters as 
close to the global minimum as possible. This is accomplished by the fast 
global convergence of DIRECT. Then, a few LM iterations are sufficient to 
rapidly refine the value of the minimum. 



generalizing this method to higher dimensions because of the difficulty in es- 
timating the Lipschitz constant. Actually, the function is often not Lipschitz 
continuous. The DIRECT algorithm overcomes these problems because it does 
not require knowledge of the Lipschitz constant. It uses all possible values of 
such a constant to determine if a region of the space of parameters should be 
broken into subregions because of its potential interest. Fig.[T]shows a schematic 
illustration of the subdivision process for a function of two parameters. 

One of the main drawbacks of deterministic global optimization methods is 
that they present a very good global convergence but a poor local convergence. 
Therefore, they quickly locate the region of the space of parameters where the 
global minimum is located, but the refinement of this solution takes a very long 
time. Consequently, the DIRECT method is a perfect candidate for its applica- 
tion as an estimator of the initial value of the parameters that are then rapidly 
refined by the LM algorithm. Since the initial point is very close to the min- 
imum, the LM method, thanks to the quadratic convergence behavior, takes a 
few iterations to rapidly converge towards the global minimum. We have con- 
firmed that this scheme reduces the number of evaluations of the merit function. 
It is important to point out that, because the most time-consuming part of the 
optimization procedure is the evaluation of the complex forward problem (with 
the solution of the statistical equilibrium equations and the radiative transfer), 
it is fundamental to minimize the number of merit function evaluations. 

A critical and fundamental point in the optimization of functions is the 
stopping criterion. In the few (unrealistic) cases where the value of the function 
at the global minimum is known (noise-free cases), it is possible to stop the 
convergence process when the relative error is smaller than a fixed value. In 
the general case with DIRECT, very good results are obtained when the hyper- 
volume (in the parameter space) where the algorithm is looking for the global 
minimum has been reduced by a given factor, / ~ 0.001, with respect to the 
original hypervolume. Another possible option is to stop after a fixed number 
of evaluations of the merit function. 

The full inversion scheme, shown schematically in Table [U is started with 
the DIRECT method to obtain a first estimation of the thermodynamical pa- 
rameters by using only Stokes /. After this initialization, some iterations of the 
LM method are carried out to refine the initial values of the thermodynamical 
parameters obtained in the previous step. Once the LM method has converged, 
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Figure 2. Example of the best theoreti cal fit obtained with the inversion 
code to the observed profiles presented in iMerenda et al.l (|2006|) . The emer- 
gent Stokes profiles have been obtained by assuming an optically thin plasma. 
The inferred magnetic field vector is given by B = 26.8 G, 6b — 25.5° and 
Xb = 161°. Furthermore, the inferred thermal velocity is Vth — 7.97 km 
s" 1 . All these value s are in close agreement with the results obtained by 
IMerenda et al.l (|2006[ K The positive direction of Stokes Q is the one parallel 
to the solar limb. 



the inferred values of u t h; " mac , and a (and At when it applies) are kept fixed. 
In a next step, the DIRECT method is used again for the initialization of the 
magnetic field vector (B, 6b and xb)- According to our experience, the first 
estimation of the magnetic field vector given by the DIRECT algorithm is typ- 
ically very close to the final solution. As a final step, some iterations of the 
LM method are performed to refine the value of the magnetic field strength, 
inclination and azimuth until reaching the global optimum. Obviously, if any 
parameter is known in advance, we can keep it fixed during the whole inversion 
process. 



4. Illustrative Applications 

Here we have only space to demonstrate the application of this new plasma 
diagnostic tool to the case of a prominence and to that of an emerging magnetic 
flux region. For more details and applications see Asensio Ramos &i Trujillo 
Bueno (in preparation). 
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Figure 3. Emergent Stokes profiles of the He I f 0830 A multiplet in an 
emerging flux re gion. The circles represent the observations shown by 
lLagg et all |2004) in their Fig. 2. The dotted lines present the best theo- 
retical fit obtained neglecting the influence of atomic level polarization, which 
corresponds to B — f070G, 8b — 86° and xb — — 160°. On the contrary, 
the solid lines show the best fit when the effect of atomic level polarization 
is included, giving B = f009G, B = 91° and X b = -161°. The best fit to 
the observations is provided by the solid lines, which indicates the presence 
of atomic level polarization in a relatively strong field region. 



4.1. A polar crown prominence 



For illustrative purposes, we have applied our inversion code to infer the mag- 
netic field vector of the pol ar crown prominence t hat generated the Stokes pro- 
files presented in Fig. 9 of iMerenda et all ( 20061 ). We have assumed that the 
prominence plasma was optically thin. Our inversion code was used to infer the 
value of the thermal velocity «th and the magnetic field vector (B, 9b, Xb)- The 
height of the promin ence atoms was fixed at h = 20", the same value used by 
Merenda et al.l (2006). After the four steps summarized in Table [IJ we ended up 
with a thermal velocity of = 7.97 km s -1 and a magnetic field vector charac- 
terized by B = 26.8 G, 6b = 25.5°, and xb = 161.0°. The inferred magne t ic fiel d 
vector is in very good agreement with that obtained by I Merenda et a~ (2006), 
namely B = 26 G, 9b = 25°, and xb = 160.5°. Therefore, we fully support their 
conclusion of nearly-vertical fields in the observed polar-crown prominence. 
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4.2. An emerging Magnetic flux region 

As pointed out bv lTruiillo Bueno Sz Asensio Ramosl ( 20071 ). the modeling of the 
emergent Stokes Q and U profiles of the He I 10830 A multiplet should be done 
by taking into account the possible presence of atomic level polarization, even 
for magnetic field strengths as large as 1000 G. An example of a spectropolari- 
metric observation of an emerging magnetic flux region is shown by the circles 
of Fig. [3] The so lid lines show the best theoretical fit to these observations of 



Lagg etal.1 (|2004f ). Here, in addition to the Zeeman effect, we took into account 
the influence of atomic level polarization. The dotted lines neglect the atomic 
level polarization that is induced by anisotropic radiation pumping in the solar 
atmosphere. Our results indicate the presence of atomic level polarization in 
a relatively strong field region (<~-d000 G). However, it may be tranquilizing to 
point out that both inversions of the observed profiles yield a similar magnetic 
field vector, in spite of the fact that the corresponding theoretical fit is much 
better for the case that includes atomic level polarization. 
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